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Abstract 

We study forward dijet production in dilute-dense hadronic collisions. By considering 
the appropriate limits, we show that both the transverse-momentum-dependent (TMD) and 
the high-energy factorization formulas can be derived from the Color Glass Condensate 
framework. Respectively, this happens when the transverse momentum imbalance of the dijet 
system, kt, is of the order of either the saturation scale, or the hard jet momenta, the former 
being always much smaller than the latter. We propose a new formula for forward dijets 
that encompasses both situations and is therefore applicable regardless of the magnitude 
of kt- That involves generalizing the TMD factorization formula for dijet production to 
the case where the incoming small-x gluon is off-shell. The derivation is performed in two 
independent ways, using either Feynman diagram techniques, or color-ordered amplitudes. 



1 Introduction 


Forward particle production observables in proton-proton (p+p) and proton-nucleus (p+A) 
collisions at the Large Hadron Collider (LHC) offer unique opportunities to study the dynamics 
of QCD at small x. and in particular the non-linear regime of parton saturation |1|. Indeed, in 
high-energy hadronic collisions, forward particle production is sensitive only to high-momentum 
partons inside one of the colliding hadrons, which therefore appears dilute. By contrast, for 
the other hadron or nucleus, it is mainly small-momentum partons, whose density is large, that 
contribute to the scattering. Such processes, in which a large-x projectile is used as a probe 
to investigate a small-x target, are sometimes called dilute-dense collisions. Since the high-x 
part of the projectile wave function is well understood in perturbative QCD, forward particle 
production is indeed ideal to investigate the small-x part of target wave function. This is true 
both in p+p and p+A collisions, although using a target nucleus does enhance the dilute-dense 
asymmetry of such collisions. 

The separation between the linear and non-linear regimes of the target wave function is 
characterized by a momentum scale Q s (x), called the saturation scale, which increases as x 
decreases. Dilute-dense collisions can be described from first principles, provided Q s 3> Aqcd- 
This condition is better realized with higher energies (as they open up the phase space towards 
lower values of x), and with nuclear targets (since, roughly, Q s ~ A 1 / 3 ). Over the years, the Color 
Glass Condensate (CGC) effective theory 2 has emerged as the best candidate to approximate 
QCD in the saturation regime, both in terms of practical applicability and of phenomenological 
success |3|. In this paper, we focus on forward dijet production in p+A and p+p collisions. 
We note that the CGC approach has been very successful in describing forward di-hadron 
production at RHIC [4j-[6], in particular it predicted the suppression of azimuthal correlations in 
d+Au collisions compared to p+p collisions [ 7 ] , which was observed later experimentally [8][9]. 

With forward dijets at the LHC however, the full complexity of the CGC machinery is 
not needed. Indeed, for the di-hadron process at RHIC energies, no particular ordering of the 
momentum scales involved is assumed in CGC calculations, while, at the LHC, the presence 
of particles with transverse momenta much larger than the saturation scale clearly must imply 
some simplifications. On the flip side, there will be other complications since further QCD 
dynamics, which is not part of the CGC framework but which is relevant at large transverse 
momenta, must also be considered. There are three important momentum scales in the forward 
dijet process: a typical transverse momentum of a hard jet, Pt, whose precise definition will be 
stated in the next section; the transverse momentum of the small-x gluons involved in the hard 
scattering, A+ and the saturation scale of the small-x target, Q s . Clearly, P t is always one of the 
hardest scales, and it is much bigger than Q s , which is always one of the softest scales. Then, 
depending on where kt sits with respect to these two, three different regimes can be defined. 

A first regime, with Q s <C kt ~ Pt, corresponds to the domain of applicability of the so- 
called high energy factorization (HEF) framework 10,11 , in which the description of forward 


dijets involves an unintegrated gluon distribution for the small-x target, along with off-shell 
hard matrix elements. That is explicitly shown in this work, starting form CGC calculations. 
While such a factorization does not occur when non-linear saturation effects are accounted for, 
we shall see that taking the Q s <C kt ~ Pt limit is tantamount to restricting the interaction 
with the small-x target to a two-gluon exchange, therefore allows to indeed write all the CGC 
correlators in terms of a single gluon distribution. Doing so, the matrix elements of the HEF 
framework are exactly recovered. 

A second regime, with kt ~ <C Ft, is where the so-called transverse momentum depen¬ 
dent (TMD) factorization [l2|| is valid. It involves on-shell matrix elements but several uninte¬ 
grated gluons distributions. In this regime, non-linear effects are present, and in the large-IV c 


2 





limit, equivalence with CGC expressions was shown in 13 . In particular, in that case the 


description of forward dijets involves only two independent unintegrated gluons distributions, 
each of which can be determined in various other processes [l4j. In the present work we shall 
keep N c finite, implying, as we show below, that a total of six independent unintegrated gluons 
distributions are needed. 

Finally, the intermediate regime Q s <C kt Pti which is naturally obtained from the two 
others by taking the appropriate limits, corresponds to the collinear regime, with on-shell matrix 
elements and the standard integrated gluon distribution. 

Separately, the HEF and TMD approaches to dijet production have been extensively studied 
in the literature 


11,15-18 and 112,19 


, but little connection has been made between them 
so far. The first result of this paper is to reveal that connection, in the context of dilute-dense 
collisions, and to show that, in fact, they are both contained in the CGC description. However, 
as already mentioned, using the CGC approach is unnecessarily complicated and one should take 
advantage of the fact that Pt 3> Q s to simplify the theoretical formulation. The second result 
of the paper is precisely to develop a new formula for forward dijets in dilute-dense collisions 
that encompasses all three situations described above, meaning that it is applicable regardless of 
the magnitude of kt- As explained below, this is obtained by extending the TMD factorization 
framework, more precisely by supplementing it with off-shell matrix elements. 

Note that the derivation of our new unified formula is performed in two independent ways: 
first using the standard Feynman diagram technique, and second by exploiting the so-called 
helicity method that employs color-ordered amplitudes 126 j. With this second method, the 
gauge invariance of the results is explicit, and the method will also prove very useful in the 
future, when processes with more particles in the final state are considered. As is the case in the 
CGC framework, our new formulation contains all the relevant limits, but it has the advantage 
that it is more amenable to phenomenological implementations than CGC calculations. In 
addition, it is also better suited to be supplemented with further QCD dynamics relevant at 
high Pt , such as Sudakov logarithms 127,28 or coherence in the QCD evolution of the gluon 


density 29-31 . These tasks are left for future work. 


The plan of the paper is as follows. In Section 2, we introduce kinematics and notations, and 
briefly present the HEF and TMD frameworks. In Section 3, we show that the HEF framework 
can be derived from CGC calculations, when the Q s <C kt ~ Pt limit is considered; namely we 
explain how the various CGC correlators reduce to a single gluon distribution in that limit, and 
show that the off-shell matrix elements of the HEF framework are indeed emerging. Section 4 
is devoted to the kt ~ Q s Pt limit, the derivation of the TMD factorization formula for 
forward dijets given in |14i] is recalled, and extended to the case of finite N c , implying six 
independent unintegrated gluons distributions instead of two. The hard factors of the TMD 
framework are computed again in Section 5, but keeping the small-a; gluon off-shell, which leads 
us to our new unified formula for forward dijets in p+A collisions. In Section 6, both the TMD 
factorization formula and the off-shell hard factors are derived again, but using color-ordered 
amplitudes, instead of Feynman diagram techniques. Finally, Section 7 is devoted to conclusions 
and outlook. 


2 Forward dijets in p+A collisions 

We shall discuss the process of inclusive dijet production in the forward region, in collisions of 
dilute and dense systems 


p(Pp) + Mpa ) -t ji(pi) + J 2 O 2 ) + X . 


( 2 . 1 ) 
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Figure 1: Inclusive dijet production in p+A collision. The blob H represents hard scattering. 
The solid lines coming out of H represent partons, which can be either quarks or gluons. 


The process is shown schematically in Fig[lj The four-momenta of the projectile and the target 
are massless and purely longitudinal. In terms of the light cone variables, iF = (v° ± v 3 )/y/2, 
they take the simple form 

Pp=]fl( l.Ot.O), PA = Jl(0,0 t ,l), (2.2) 

where s is the squared center of mass energy of the p+A system. 

The energy (or longitudinal momenta) fractions of the incoming parton (either a quark or 
gluon) from the projectile, aq, and the gluon from the target, aq, can be expressed in terms of 
the rapidities and transverse momenta of the produced jets as 

xi = — \ P ' 2 = ~^= (\pit\e yi + \p 2 t\e V2 ) , 

Pp V s 

X 2 = Pl + _ P2 = - 7 = {\pit\e~ yi + \p 2 t\e~ V2 ) , 

where pit, P 2 t are transverse Euclidean two-vectors. By looking at jets produced in the forward 
direction, we effectively select those fractions to be aq ~ 1 and aq 1. Since the target A is 
probed at low aq, the dominant contributions come from the subprocesses in which the incoming 
parton on the target side is a gluon 

qg -+ qg, gg^ qq, gg -+ gg- (2-4) 

In dilute-dense collisions, the large-a; partons of the dilute projectile are described in terms of the 
usual parton distribution functions of collinear factorization f a / p , with a scale dependence given 
by DGLAP evolution equations. By contrast, the small-a; gluons of the dense target nucleus 
are described by a transverse-momentum-dependent distribution, which evolve towards small x 
according to non-linear equations. Moreover, the momentum k of the incoming gluon from the 
target, besides the longitudinal component k~ = X 2 y/s/ 2 , has in general a non-zero transverse 
component, kx, which leads to imbalance of transverse momentum of the produced jets 

\kt\ 2 = \pit + P 2 t\ 2 = \pit\ 2 + \p 2 t | 2 + 2\Pit\\P2t\ cos Acf >, (2.5) 


(2.3a) 

(2.3b) 
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with k^ = — \kt\ 2 . Here, by k t we mean a four-vector, as opposed to kt = Pit + P 2 t, which is 
a two-dimensional vector in the transverse plane. They are simply related by: kr = (0,fct,0). 
Using the notation defined above, the gluon’s four-momentum can be also parametrized as 

k = x 2 pA + kr ■ (2.6) 

The Mandelstam variables at the partonic level are defined as 


s= (p + k) 2 = (pi +p 2 ) 2 = 

\Pt\ 2 

z(l-z) 

t = (P2 ~ pf = (pi - kf = 

1 P2t\ 2 

1-z’ 

u = (pi- pf = (P2 ~ k) 2 = 

M 2 


(2.7a) 

(2.7b) 

(2.7c) 


with 

z = + Pl , and P t = (1 - z)p u - zp 2t • (2.8) 

Pi + P 2 

They sum up to s + t + u = kj,. 

Note that we always neglect the transverse momentum of the high-x partons compared with 
that of the low-x parton \kt\. This is justihed in view of the asymmetry of the problem, xi ~ 1 
and X 2 < 1 , which implies that gluons form the target have a much bigger average transverse 
momentum (of the order of Q s ) compared to that of the large x partons from the projectile 
(which of the order of Aqqd)- And even when the transverse momentum imbalance of the dijet 
system is of the same order as the jet transverse momenta themselves, implying that both parton 
distributions are probed in their radiative tail, the small x 2 (BFKL) evolution implies a 1 /kt 
behavior on the target side, while DGLAP evolution implies a 1 /k 2 behavior on the projectile 
side. 

To take into account small-x effects in dijet production, an approach that has been broadly 
used in phenomenological studies involves the so-called high energy factorization (HEF) for¬ 
mula [l5 


dQ-pA—t dijets+Jf 

dy i dy 2 d?p 1 1 d 2 p 2t 


16tt3(xix 2 s) 2 Xifa/ P {xi : p 2 ) \M ag *^ cd \ 2 P g/A (x 2 , kt) 1 + 5cd - 

a,c,d 


(2.9) 


This formula makes use of the unintegrated gluon distribution P g /A that is involved in the cal¬ 
culation of the deep inelastic structure functions. It is determined from fits to DIS data, and 
then used in Eq. (2.9), along with matrix elements that depend on the transverse momentum 
imbalance (2.5). Even though the high energy factorization is not strictly valid for dijet pro¬ 
duction, there exists a kinematic window, the dilute limit Q s <C \pit\, \p 2 t \1 \kt\, in which it can 
be motivated from the CGC approach. We shall demonstrate this explicitly for all channels in 
the next section. 

A second approach, valid in the regime where the transverse momentum imbalance between 
the outgoing particles, Eq. (2.5), is much smaller than their individual transverse momenta, is 
the so-called transverse momentum dependent (TMD) factorization. This limit, | pit + P 2 t\ 
\pit.\, \p 2 t\, or life* | <S \Pt\, corresponds to the situation of nearly back-to-back dijets. Even though, 
in general, there exists no TMD factorization theorem for jet production in hadron-hadron 
collisions, such a factorization can be established in the asymmetric “dilute-dense” situation 
considered here, where only one of the colliding hadrons is described by a transverse momentum 
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dependent gluon distribution. Again, selecting dijet systems produced in the forward direction 
implies x\ ~ 1 and X 2 <C 1, which in turn allows us to make that assumption. The TMD 
factorization formula reads (so far, this has been obtained in the large- N c approximation, but 
this restriction will be lifted in the present work) 1131 


at 


^Q-pA-^ dijets+X 


dy 1 dy 2 d 2 p lt d 2 p 2 t (xi x 2 s) 2 


^xi/ a/p (x!,/x 2 )^ 

CLQ—>cd Fag{x 2,b) 


1 


a,c,d 


1 + d c d 


( 2 . 10 ) 


U) 

where several unintegrated gluon distributions Tag with different operator definition are involvec 


(z) 

and accompanied by different hard factors H y a g^cd m Those hard factors were calculated in 
as if the small-X 2 gluon was on-shell (i.e. |A:*| = 0). The kt dependence survived only in 
gluon distributions. 


13 

the 


By restoring the kt dependence of the hard factors inside formula (2.10), we can make the 


bridge between the HEF and TMD frameworks and obtain a unified formulation which encom¬ 
passes both the dilute and the nearly back-to-back limit. Note that we follow the conventions 
used in earlier papers that dealt with these formalisms, such as Ref. 

|2 


15 and 13 respectively. 


the hard factors H 


(0 


of the 


Therefore, contrary to the HEF matrix elements \M. a g*^cd\ , ucm i^i ±± ag ^cd 
TMD factorization are defined without the g' 1 factor. In addition, the definition of the gluon 
distribution also differ by a factor it. The integrated gluon distribution X2 f g /A is obtained from 

f dk 2 J-g/A i n the HEF formalism, and from f d 2 kt Tag in the TMD formalism. 

Finally, let us point out that, in the frameworks described above, one emits radiation in the 
transverse direction that one has no control over, as it is part of the small-x gluon distributions 
and therefore is treated fully inclusively. To be more specific, at this level, transverse momentum 
conservation is obtained either by several particles of average transverse momentum Q s , or by a 
third hard jet, depending on the magnitude of \kt\. Due to the small-x evolution, that radiation 
is ordered in rapidity, therefore it does not contribute to the measured forward dijets systems. 


3 High energy factorization derived from CGC: 
the |p lt |, \p 2t \, |fc t | > Q s limit 

We shall demonstrate that the high-energy factorization formula for double-inclusive particle 
production, Eq. (2.9), is identical to a result obtained from the CGC formalism in the dilute 


target approximation. This is a limit where all the momenta involved in the process are much 
larger than the saturation scale: \pit\, \p 2 t\, \kt\ t$> Q s ■ Here, we show explicitly the equivalence 
of the HEF and CGC formulas for the qg* —> qg channel and only provide the final results for the 
two other channels, as the derivations proceed identically for all of them. We derive the CGC 
cross sections for the qg* —>• qg and gg* —> qq channels in the dilute limit following a procedure 
developed in Ref. j32j where only the gg* —> gg sub-process was considered. 

The amplitude for quark-gluon production is schematically presented in Fig. [2] as in Ref. [7j. 
In the left diagram, the emission of the gluon from the quark happens before the interaction 
with the target, and in the right diagram the emission occurs after the quark has interacted 
with the target. There is a relative minus sign between the two cases as explained in details 
in Ref. [7|. Multigluon interactions of quarks and gluons with a target, in the CGC formalism, 
enter as Wilson lines in the expression for the amplitude. A quark propagator is represented as 
a fundamental Wilson line, while a gluon propagator as an adjoint Wilson line. As a result, the 
cross section involves multipoint correlators of Wilson lines. In particular, the amplitude from 
Fig. [ 2 J after squaring, has four terms: a correlator of four Wilson lines, S^\ corresponding to 
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Figure 2: Amplitude for quark-gluon production in the CGC formalism. Left: the gluon is 
radiated before the interaction with the target. Right: the gluon is radiated after the interaction 
with the target. The two terms have a relative minus sign. 


interactions happening after the emission of the gluon, both in the amplitude and the complex 
conjugate, then a correlator of two Wilson lines, representing the case when interactions 
with the target take place before the radiation of the gluon in both amplitude and complex 
conjugate, and two correlators of three Wilson lines, S^ 3 \ for the interference terms. In all the 
cases the splitting function is the same, and is given by the product of the quark wave functions: 
^piPiPl ■> x/ — x — b). The total expression for the inclusive cross section in CGC 

is then given by the following formula [7]: 


da(pA—¥ qgX) = _ Z ) p + Xl f f ( Xl , p 2 ) \M{p,pi,p 2 )Y 

dy\dy2d z j)\td z p2t 


(3.1) 


where the amplitude squared, \M.(p,pi,p- 2 )\ , has the form: 
\M{p,pi,p 2 )\ 2 = 


d 2 x d 2 x' d 2 b d 2 b' 


g —jpir(x—x') e —ip 2 f(b—b') 


(27t) 2 (27r) 2 (27r) 2 (27r) 2 

X J2 $$3 & Pi ’’ x/ “ (LpI.X- b) 


Xa/3 


X {5 , g3L[b,x,b , ,x / ;x 2 ] - 5^[b,x,b' + z{x! - b');x 2 ] 

-S^[b + z(x - b),x , ,b / ;x 2 ] + S^f[b + z(x - b),b' + z(x' - b');x 2 ]| , (3.2) 


ex- 


where are mixed-space quark wave functions and are correlators of Wilson lines 
plained in details below. Following the notation from Fig. [I] and Eq. (2.8), we use the fraction 
of the plus components of four-momenta, z, with p\ being the four-momentum of the outgoing 
gluon and p 2 , the four-momentum of the outgoing quark. 

The fundamental, U(x), and adjoint, V(x), Wilson lines are defined as path-ordered expo¬ 
nentials of the gauge field (written here in the A + = 0 gauge): 


U(x) = V exp 


W 


dx + A a (x + ,x)t a 


and V (x) = V exp 


W 


dx + A a ( x + ,x)T a 


(3.3) 


where t a and T a are the generators of the fundamental and adjoint representations of SU(N) 
respectively. The traces of products of Wilson lines appearing in the cross section are defined 
in the following way: 


o( 4) 
9393 


(b,x,b', 


x) 


C f N c 


Tr 


( 'u(h)U\h')t d t c ) V (x)F t (x / 


1 cd 


X2 


(3.4) 
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(3.5) 


Sg(b,x,z') = ^ (TV (u\z')t c U(b)t d ) F cd (x)^ ; 

sgW) = ^ ( Tr (l/( z )^ t (z , ))) iB2 • (3-6) 

The CGC average is taken over the background filed evaluated at Y = ln(l/x 2 ). The product 
of wave functions in the massless limit is: 


5Z (P>Pi"> u 0^(P, Pi". u ) 

A a/3 


8 vr 2 


u • u 


Pi 


r(l + (l-^) 2 )- 


(3.7) 


Introducing a change of variables, u = x — b and v = zx + (1 — ^)b (and similar for the primed 
coordinates), we get (7j: 

\M(p,p!,p 2 )\ 2 = j (^^^2 elP4 ' (U,_U) IZ^(P,P^, u/ )A(P,Pmu) 


r d 2 v dV 

(27r) 2 (2 tt) 2 


0 ifcf(v'-v) 


'S'ogL^jX, b^x 7 ) - 5^(b,x,v') - Si^vy.b'J + S^^y) .(3.8) 


? (3) 

539 


(3 L 

993 V 


? ( 2 ), 

m 


The conjugate momentum to u' - u is Pi = (1 — z)pu — zp 2 t, and the one corresponding to 
v 7 — v is the total transverse momentum of the produced particles kt = pit + P 2 t ■ In terms of 
fundamental Wilson lines only: 

^SL(b,x, b',x') = 2 ^ (Tr (p(b)pt( b ')P(x')P t (x)) TV (V(x)pt( x ')) ( 3 . 9 ) 

--^Tr(p(b)pt( b '))\ , 

iVc ' / X2 

and 

5g(b,x,v') = <^Tr (P(b)pt (x) ) T r (u(x)rf(v')) ~ ^TV ([/(b)pt( v ')^ . 

" 2(3 ' 10) 

In the dilute target limit we allow for only up to two gluon exchanges between the Wilson 
line propagators and the nucleus. Accordingly, we expand the Wilson lines to second order in 
the background field: 

dx + A~(x + , x) — J dx + dy + V {A _ (x + ,x)A _ (y + ,x)} + 0(A 3 ) . (3.11) 


P(x) « 1 +ig 


To this order, the expectation values of the four- and three-point correlators are simply expressed 
in terms of the dipole operator SA (v, v'). The dilute target approximation gives only a leading 

result in |v — v'| 2 Q 2 for the expectation value of (v. v r ). which is equivalent to taking the 
limit \kt\ S> Q s - Similarly, when all the momenta involved in the process are much larger than 
the saturation scale, the correlators entering the cross section get the following expressions: 

S q%g( h ’X’ h ' , x ') = 1 - g 2 N c r x2 (x-x') — g 2 N< 2 N 1 r x 2 {h - b 0 

[-Tra( x - b ) + r ^( x ' -b') -Px a (x' -b) -r x2 (x-b')] ; (3.12) 
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:( 3 ), 


g 2 N c 


r Nn a 2 

2 c 7(x-v') +2 9 w 7( b-v'); 

(3.13) 

|jVc 2 _ 1 r r2 (v v 7 ). 

2 N c 2V J 

(3.14) 

< X2 (x + ,0) ~'y X2 (x + ,r)] , 

(3.15) 


qgq^ ' ’ > 2 2V ' 

s^Ow') = 1 - g‘ 

In the above equations: 

r x - 2 ( x - y) = J dx + [- 

where r = x — y and 7 a; 2 (:r , r) is related to the expectation value of the two-field correlator: 

(A~(x + ,x.)A~(y + , y)) x2 = 5 ab 8(x + - y + )'y X2 (x + , x - y) . (3.16) 

Using the expressions for the multi-point functions SW, we get the following result for the 
amplitude squared: 


\M(p,pi,p 2 )\ =4tt g N c (l + (1 - z) W 


d 2 u dV ei p t . (u /_ u) u • u 7 


J ( 27 r ) 2 ( 27 r ) 


u| 2 |u 7 | 2 


X / v) [r x2 (x - b 7 ) + U X2 (x - b) + r x2 (x - v) 


( 27 t ) 2 ( 27 t ) 2 


j\r 2 - 1 

+r X2 (v - x) - 2r x . 2 (x - x) —r x . 2 (b - b') 
AT 2 — 1 1 1 

~iv 2 -^(v- v) - ^r X 2 (b- V) - ^r X 2 (v-b') 


(3.17) 


We perform the integrals in the above expression by changing the variables from v and v 7 to r 
and B. The integrals over the transverse distances of the type r = v — v 7 are equivalent to the 


Fourier transform of Eq. (3.15) and give the unintegrated gluon distribution: 


fx 2 (h) = -kt / d 2 r r x2 (v)e lkt ' r = kf / dx + 'y X2 (x + , k t ) . 


(3.18) 


In our approximation, the correlators do not depend on the impact parameter B = (v + v 7 )/ 2 . 
The integrals over B factorize and give the transverse area of the target: j d? B = S±. Finally, 
the rest two integrations reduce to: 


d 2 ue~ iPtU 


u 


= —27ri¬ 


ll 


Pt| 2 ' 


(3.19) 


In terms of the unintegrated gluon distribution, the amplitude squared then gets the form: 


\M(p,pi,p 2 )\ = 


(2tt)‘ 




fx 2 (k t 

kt 


(1 +( 1 - 7)4 

Pi 


(■ N 2 ~ 1 ) 1 (■ N 2 - 1) 1 1 1 Pi ~ Pit Pi • P2t Pit • P2t 

2 IV 2 P 2 2IV 2 p\ t pi IV 2 P 2 pf t P}p\ t p\ tV \ t J 


, (3.20) 


We want to show that Eq. (3.20) reproduces the HEF formula (|2.9[) with the appropriate 


unintegrated parton distribution function and off-shell matrix elements. For this purpose, we 
need to find a relation between the unintegrated gluon distribution used in the above equation, 


f X2 (kt ), and P g /A(x 2 ,kt), which appears in the HEF formula (2.9). This is easily done by 
considering the deep inelastic scattering process, since P 9 /a( x 2 , kt) is precisely the unintegrated 


9 

























gluon distribution involved in the formulation of the 7 * + A —> X total cross section, and is 
therefore related to the qq dipole scattering amplitude in a straightforward manner (see for 
instance 116,331): 


N r . 


FglA(x2,kt) as ( 2 7T)3 


d 2 vdV e -iM v - v ') V 2 


1 - sff( v > v/ ) 


In the weak-field limit, using formula (3.14), this gives the relation 

47T 2 


f(kt) Sj_(N 2 _ -Q Pg/A(x2, kt) 


(3.21) 


(3.22) 


Then, the cross section for the qg production channel from Eq. (3.1) can be written in a more 
compact form 


dyidy2d 2 Pud 2 L 2 F Xlfq ^ z)P ^ z) 


1 + 


(l-z) 2 Pl 2 t 1 z 2 p 2 2 t 


p? 


N 2 if j 


Fg/A(x2,h) 


Pit P2t 


where Pgq(z) is related to the quark-to-gluon splitting function and is given by: 

• 


(3.23) 

(3.24) 


It turns out that the above expression for the quark-gluon production cross section is identical 

-2 

to the result in the HEF formalism, Eq. (2.9), containing the off-shell amplitudes \M a g*^cd\ ■ 
The latter have been calculated in Refs. (Uf [34] and (35]. 

The equivalence of the CGC and HEF formulas in the dilute limit can be shown in a similar 
way for the cross sections of the other two subprocesses, gg* -A qq and gg* -A- gg. The CGC 
results for the cross sections in this limit are: 

_ J_ (l~z ) :2 Pit +Z 2 p 2 t1 Pg/A{x 2, kt) 

m p? 


da(pA -A qqX) 


at 


dyidy2d 2 pitd 2 p2t ^C F ir 


Xlfg/p{xi, p 2 )z(l-z)P qg (z) 


Pit P^. 


and 132 

dcr(pA -A ggX) 


asNC Xlfg/ p {Xl, p 2 )z(l-Z)Pgg{z) 


dyidy2d 2 p u d 2 p2t nC F 
The expressions for Pqgi z ) and Pgg(z) have the form: 

Pqg(z) = z 2 + (l-z) 2 , Pgg(z) = 


1 + 


(3.25) 

(l-z) 2 pf + z 2 p£ 1 Xg/ A (x 2 , h) 


1 - 2 


Pi 2 


1 — z , 

+- +z(l-z). 


pit pi 


(3.26) 

(3.27) 



in order to investigate the so-called geometric scaling regime, where saturation effects are felt, 
even though Q s <C kt- The full saturation region, Q s ~ kt, is however, in principle, out of reach 


of formula (2.9). Along these lines, an estimate of saturation effects was obtained in [42.43 


using the BK equation extended to include the same higher-order corrections as included in the 


linear case 39 . 
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4 TMD factorization for nearly back-to-back jets: 

the |p lt |, |p 2t | > \k t \,Q 8 limit 


In this section we discuss the special case of nearly back-to-back jets, | pu + P 2 t\ \pit\i \P 2 t\, 


. . . . . (i) 

where the differential cross section is given by formula (2.10). Several gluon distributions Fag, 
with different operator definition, are involved here. Indeed, as explained in [l2|, a generic 
unintegrated gluon distribution of the form 

F(x 2 ,k t ) n ^ e 2 [ ^ )Xr { F i- f ^ F i- (0) ] \A) , (4.1) 

J (27 t) 6 p a 


where F l ~ are components of the gluon field strength tensor, must be also supplemented with 
gauge links, in order to render such a bi-local product of field operators gauge invariant. 

The gauge links are path-ordered exponentials, with the integration path being fixed by the 
hard part of the process under consideration. Therefore, unintegrated gluon distributions are 
process-dependent. 

In the following, we shall encounter two gauge links U M and IA^~\ as well as the loop 
U PI = if WifHt = IjF'U +]t. These links are composed of Wilson lines, their simplest expression 
is obtained in the A + = 0 gauge: 


= 17(0, Too; 0)Z7(±oo, £ + ; £) with U(a, 6; x) = V exp 


ig / dx + A a (x + , x)f a 


(4.2) 


but the expressions of the various gluon distributions given below are gauge-invariant. From now 
on, F l ~ (£ + ,£) is simply denoted as F(£), and the hadronic matrix elements (A|...|A) —> (...). 
Note however that they are different from the CGC averages (■ ■ ■) of the previous section. 
Indeed, the normalization of the hadronic state |A) is defined as (A'\A) = (2-k) 3 2p~^ 5(p^ — 
Pa) ijP At ~ p'At)i w hile the CGC averages are normalized as (1) = 1. As explained in 1131, 

the two can be related by making the replacement (• ■ - ) X2 

This approach to dijet production in proton-nucleus collisions was analyzed in Ref. 
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The TMD factorization formula ( |2. 10 ) was derived there in the large- N c limit, and shown to 
be equivalent to CGC calculations (e.g. formulas ( |3.1[ ) and (3.2) in the case of the qA —> qg 
channel), after taking the limit |pit|,|p2t| ^ \kt\,Qs- In this section, we derive the TMD 
factorization formula keeping N c finite. We obtain corrections to the hard factors H^^ cd 
previously derived, and we calculate new hard factors corresponding to gluon distributions that 
were omitted before (as they were vanishing in the large- N c limit). The finite N c extension 
prevents one to make a further simplification, called correlator factorization, essential to relate 
the TMD factorization and the CGC formalism, but gives completeness to the main result of this 
paper, i.e. the new factorization formula we propose below is valid for finite N c . We also check 
explicitly the gauge invariance of these hard factors by computing them in a gauge different 
from the one used in 
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An important fact to note is that, as a consequence of the | A:* | <C |pit|, |f> 2 t| limit, the kt 


dependence in (2.10) survives only in the gluon distributions, and the hard factors are calculated 


as if the small-X 2 gluon was on-shell. That is, looking at the hard partonic interaction represented 
by the blob H in Fig. [lj /c 2 = — \kt\ 2 is set to zero, and s + t + u = 0. 


4.1 The qg —> qg channel 

The complete set of independent cut diagrams contributing to this channel is shown in Fig. [3] 
(mirror images of diagrams (3), (5) and (6) give identical expressions). 
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Figure 3: Diagrams for qg 
identical contributions. 


qg subprocess. The mirror diagrams of (3), (5) and (6) give 


The cross section for a quark-gluon scattering involves only two different TMD gluon distri¬ 
butions as given in Ref. 


12 


d(J pA^qgX 


at 


d 2 P t d 2 k t dy 1 dy 2 {x\x 2 s) 2 


xif q/p ( Xl 


i =1 


(4.3) 


with: 


7r(l) 
J 99 


77(2) 
J 99 


= 2 


dt + d 2 a 

(2tt ) 3 p~ a 


(2tt) 3 Pa 


-e iX2 p A t+-ikft /-jy \f^U^F( 0)U^}) = x 2 G {2 \x 2 ,k t ) , (4.4) 

a ' L J / 

(4.5) 


2 f d£ + d £ c ix2P A Z + -ikft / rjy 


F (0 Tl ^ ] W +]t F (0) ZY [+1 


iV r 


These are the same gluon distributions as in the large- N c limit 1131, no additional ones are 


present in this channel. The only difference in the expression (4.3) when we go to finite N c will 


appear in the hard factor Hqg\ qg associated with Fqg . That gluon distribution is sometimes 
also denoted x 2 G^ 2 \ and is called the dipole distribution , since it is the one that enters the 
formulation of the inclusive and semi-inclusive DIS. 

In the CGC approach, x 2 G^ can be related to the qq dipole scattering amplitude, and there¬ 
fore linked to the gluon distribution used in the HEF formalism: P g /A{x 2 , kt) = irx 2 G( 2 ' 1 (x 2 , kt). 
That distribution is not sufficient however to compute the forward dijet cross section when 
\kt\ ~ Q s (*.e. the case considered in this section). For completeness, we note that a detailed 
derivation of this relation between formula ( 3.21| ), involving a CGC correlation function, and 
formula (4.4), involving matrix elements defining TMDs, can be found in Appendix A of 


13 


The exact results for the two hard factors read 


HqgUqg ~ ~ N 2 - 1 ^ + 2j ° 5 + 2j ° 6 ’ 

H £l q9 = ^i + ^^ + 2D 3 , 


(4.6) 

(4.7) 


where DiS are the squared and interference diagrams corresponding to the qg -A qg channel, 
following the numbering of Fig. [3j Each term Di = C Ui h{ represents the product of the color 
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(i) 

factor, C Uj , and the hard coefficient, hi. What kind of diagrams enter the hard factors H qg ^, qg 
depends on the type of the gauge links appearing in each of them. As summarized in table IV 
of Ref. 12 , the distribution J- q ' g is present in diagrams (1), (2), (4), (5) and (6), while the 
distribution T qg appears in diagrams (1), (2) and (3). The Di components were computed 
in Ref. 1131 (table II) in an axial gauge with the axial vector, n, set to n = p, for both the 


incoming and the outgoing gluon, where p is the four-momentum of the incoming quark, as 
defined in Fig. [l] Formulated differently, the polarization vector of each external gluon was 
chosen such that, besides with the momentum of the gluon, their inner product with p vanishes. 
We recovered the same results for DjS in that gauge and performed the same calculation in a 
different gauge with the axial vector set to n = p for the incoming gluon and n = p 2 for the 
outgoing gluon [^] The results for the hard factors H q l\ qg and H q ^l^ qg at finite N c are identical 
in both gauges and they read 


qg^qg 


#( 2 ) 

qg^-qg 


u(s 2 + u 2 ) 1 (s 2 + u 2 ) 

2 st 2 + 2 N 2 su 

s(s 2 + u 2 ) 

2 iii 2 


(4.8) 


(4.9) 


The hard factors and the TMDs entering the factorization formula (4.3) are all gauge invari¬ 


ant. In principle, that leaves us some freedom and the factorization formula can be rewritten 
with new hard factors and the corresponding new gluon distributions formed as linear combi¬ 
nations of the the old ones. 

For reasons that shall be discussed in detail in Section [6j let us define the new hard factors 
for the qg —> qg subprocess 


iff 1 ) = iff 1 ) + — #( 2 ) 
qg^qg qg^qg ~ jy2 qg^-qg 

and the corresponding new gluon TMDs 


and iff 2 ) = ^_ -H f 2 ) 

d,nu iy q g —>qg N 2 11 qg^qg ■ 


(4.10) 


qg^qg 

— J 7 ) 1 ) 

J qg > 


(4.11) 

<p( 2 ) 

qg^qg 

_ 1 t 

' f) _)_ ]\[ 2 jz(2)l 

J qg ~ qg ) > 

(4.12) 

N 2 - 1 1 


such that the factorization formula (4.3) now takes the form 
da'PA^qgX 


( Xl x 2 s) 2 


iff 1 ) <p( 2 ) iff 2 ) 

qg^qg qg^qg ^ qg^qg qg^qg 


d 2 P t d 2 k t dy l dy 2 

The explicit expressions for K qg ’^ qg and K qg ^ qg are given in Table [1 

4.2 The gg —> qq channel 

The independent cut diagrams contributing to this channel are shown in Fig. [4l 
In addition to the two gluon distributions, Pgl/ and P gg \ used in Ref. 1131 


(4.13) 


the result to 


12 


44 


J-gg (also sometimes denoted X 2 and 


all orders in N c involves a third distribution 
called the Weizsacker-Williams gluon distribution ), and the differential cross section reads 

da pA ^w x 9 3 


at 


d 2 P t d 2 k t dyidy 2 (aq x 2 s) 2 




gg^-qq 


(4.14) 


1=1 


J The choice of axial gauge vectors for external gluons corresponds to the choice of the reference momentum 
for their polarization vectors, see for example [26], and is arbitrary for gauge invariant quantities. Thus, the 
independence on those gauge vectors can be used to confirm that the result is gauge invariant. 
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K ( x ) K^ 

ag^cd ag^f-cd 


qg qg 

gg qq 


s 2 + u 2 l 2 s 2 — t 2 ' 

- fr A - 

2 t 2 su i N 2 


Cp sis 2 + u 2 ) 
N c t 2 u 


1 (t 2 + u 2 ) 2 1 i 2 + u 2 

2 N c § 2 tu ~ 2 C f N 2 s 2 


gg ->• gg 


2N C (s 2 — tu) 2 (P + u 2 ) 
Cf t 2 u 2 s 2 


2iV c (s 2 — iti) 2 
Cf tus 2 


Table 1: The “new” hard factors following from simplified effective TMD factorization of 
Eqs. (4.13), (4.25) and (4.52) in the case with all partons being on shell. 


XL/ 


/ i \ 

(i) 




Figure 4: Diagrams for gg —> qq subprocess. The mirror diagrams of (3), (5) and (6) give 
identical contributions. 


with the three gluon TMDs defined as 

= 2 [ f * e iX2P ~^ + ~ ikt < ( Tr 


(2vr ) 3 p A 

J&) = 

99 ' ( 2ir) 3 p A 


(0 Tr ^ [ % Wf(0)uM 

f d ^ + f{ e iX2P ^ + ~ ik ^A- (t. r \f (fl^tl Tr |f(0)Z7 [ /\ 
J {2ttYp a N c \ L J L J / 


(4.15) 

(4.16) 


Fp) = 2 


dtd 2 d 


99 ~ 1 (2 n) 3 p A 


e ix 2 P A Z+-ik t -s ^ty F(O^ [+lt F(0)^W ^ =® 2 G( 1 >(® 2 ,fc t ) . (4.17) 


The appropriate hard factors are constructed from the expressions corresponding to the diagrams 
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(l)-(6) depicted in Fig. [4j using the following formulas 


H, 


(i) 


H, 


gg^qq 

( 2 ) 


-D\ + „ D'2 + + 2D§ + 2Z?6 


2Cf 


2C f 


II, 


99^99 

(3) 


— —2iV c -D 3 — D 4 — 2 D 5 — 2F>6 , 


1 


99^99 


N? - 1 


Di - 


1 


TV? - 1 


D 2 + 2 D 3 . 


(4.18) 

(4.19) 

(4.20) 


Again, the components Di = C Ui hi were computed in 1131 (table III) and they were used there to 

( 12 ). . 

determine the hard factors Hg g ’A qq in the large N c limit. Here, we generalize the results of 


13 


to the full, fmite-lV c case. The calculation can be most readily done by exploiting crossing 
symmetry that relates the qg qg and gg —> qq channels. This allows for identification of 
the diagrams between Figs. [3] and [4] and enables one to recycle the Di expressions calculated 
in the previous subsection. For example, the expression corresponding to the diagram (1) from 
Fig-i with the incoming and the outgoing legs connected, is identical to the already computed 
expression for the diagram (4) from Fig. [ 3 ] (modulo a color averaging factor and swapping of 
the momenta p\ «->• p). Similarly for all the other diagrams. That gives the following set of hard 
factors for the gg —> qq subprocess: 


H, 


(i) 


1 ( t 2 + u 2 ) 


2\2 


II, 


99^99 

( 2 ) 


II, 


99^99 

(3) 


4 C F 

s 2 ut 

Ek 

i-H 

+ u 2 

2Cf 

s 2 ’ 

1 

P + u 2 


99^99 


4 N 2 C F tu 


(4.21) 

(4.22) 

(4.23) 


Of the three hard factors, Hg g ^ qq , only two are independent. The third hard factor, H gg \ qq , 
can be expressed a^] 

( 3 ) _ 


H, 


99^99 


1 (tt( 1) , rr(2) 

jy 2 y^gg^qq “T 11 gg^fqq 


(4.24) 


Therefore, the cross section for quark-antiquark production can be rewritten with only two hard 
factors and two gluon distributions that are linear combinations of F q ^, and F q 'g : 




da pA ^ gqX 

d 2 P t d 2 k t dyidy 2 (aqx 2 s ) 2 Xi ^/p^ X i ^ ^ 


-K^\ - + ^ (2 L - 

gg^qq gg^qq ~ gg^qq 99^qq 


In the above, we defined the new gluon TMDs as 

■if 1 ) — _ I _ ( at 2 tt(i) _ 7t(3A 


gg^qq 

( 2 ) 

99^99 


]y2 _ l \N C Fgg - J gg J , 

= —]\[ 2 J^( 2 ) 4 . jr{ 3) 

■ L1 c u gg 1 u a a ? 


99 


-(*) 


and the hard factors Kg q _^ qq as: 

jA 1 ) _ — ~~ ^ tt{ 1) _ j 1^(2) __ 3 _/r( 2 ) 

99^qq at 2 gg^-qq gg^-qq i\r2 gg^-qq • 

Jv c c 

The explicit expressions for the latter are given in Table [l] 


(4.25) 

(4.26) 

(4.27) 

(4.28) 
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Figure 5: Set of diagrams for the gg —> gg subprocess involving only 3-gluon vertices. The 
mirror diagrams of (3), (5) and ( 6 ) give identical contributions. 


4.3 The gg —> gg channel 

Finally, the independent cut diagrams for the gg —> gg channel are given in Figs. [5] and [ 6 j and 
the corresponding differential cross section for two-gluon production reads: 


da pA^ggX 

d 2 P t d 2 k t dyidy 2 


(® 1 ® 2 S ) 2 Xlf 9/p(' Xl ^ 


2\ f(i) u(i) 
> / j J gg gg- 


2—1 


99 gg->gg 


(4.29) 


(12 3). 

The J-gn ’ ; distributions are the same as the ones introduced in the previous section in Eqs. (4.15 )- 


(4.17). The remaining three are 12 


•4? = 2 


*0? = 2 


= 2 


d £ iJ.ii f + —jA'. f 

( 27 t) 3 Pa 

d£ d~£ c iX2P~^(. + -ikf$. 

(2vr fpl 

d£ d £ Jx 2 p^ + -ikf^ 


TV 


TV 


TV 


F(^)U [ -^F(Q)U [ - 
F (£) U^U [+ ^F (0) WPIW W 


{2tt) 3 Pa 

The associated hard factors are constructed a^J 


F (0 ^ [+lt F (0)£7 [+1 


Tr [2/P1] Tr [wR] \ 

““iVc A^T” / 


i/W 

gg^gg 

H ( 2 ) 

gg^gg 


— 2^ >1 + 2 ’ 
= 2H3 — -D4 — 2H5 — 2D 6 , 

r 2 

b - ,. „ 

gg^gg c gg->gg 




(4.30) 

(4.31) 

(4.32) 

(4.33) 

(4.34) 

(4.35) 


The calculation of the gg —>• gg subprocess requires inclusion of diagrams with four-gluon 
vertex. Therefore, in general, the expressions D. t in the above equations contain contributions 


The same relation holds of course already at the level of Eqs. (4.18 l-l 4.201. 

(6) 


3 Note that what is called H gg ^ gg in Ref. 


and H. 


(6) 


13 


is now H), 


. Out of six hard factors, only H^ g ^ gg , H gg ^ gg 


r( 2 ) 


survive in the large- N c limit. 
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(7) 





Figure 6: Set of diagrams for the gg —> gg subprocess involving 4-gluon vertex contributions. 
The mirror diagrams of (8), (9) and (10) give identical contributions. 



Figure 7: Color indices for the cut four-gluon squared matrix element. 


from both, the 3-gluon and 4-gluon vertex diagrams, the latter shown in Fig. [6j The correspond¬ 
ing expressions were computed in |13|, where they were used to determine the hard factors in 
the large- 1V C limit. Below, we generalize the result of Ref. 13 to the case of finit e-N c , with the 
help of the exact definitions given in Eqs. ( 4.33[)- (4.35). The six hard factors read 


i7 (1) 

a 9^ a a 

#( 2 ) 

99^99 


#( 6 ) 

99^-99 


N c ( t 2 + ii 2 )(s 2 — tu) 2 
Cf u 2 t 2 s 2 

2N C (s 2 — tu) 2 
Cf uts 2 



rC) 

99^99 


N 2 H (4) 

c ±± gg^gg 


N 2 H (5) 

c 99^-99 


N c (s 2 — tu) 2 
Cf u 2 t 2 


(4.36) 

(4.37) 

(4.38) 


To get further insight into the above results, we have performed an independent calculation 
in a gauge with non-vanishing 4-gluon vertex contribution, with the axial vectors defined as: 


n = p for the gluon k , n = k for the gluon p , 

n = P 2 for the gluon p \, n = p\ for the gluon p 2 ■ 


(4.39) 


The contributions to Di s in this gauge, coming from diagrams with 3-gluon vertices only and 
depicted in Fig. [5j are given in Table [2] 


17 
















ft! 3 ’ 

Q 

(1) 

4s 6 + 4 ts 5 + 17 Ps A + 36t 3 s 3 + 24 t A s 2 + 8 t 5 s + 4 1 6 

N c 

s 4 i 2 

2 C F 

(2) 

s 6 + 2 ts 5 + 33 f 2 s 4 + 60t 3 s 3 + 44t 4 s 2 + 16t 5 s + 4t 6 

N c 

s 4 (s + 1) 2 

2 C F 

(3) 

2s 6 — 9 ts 5 + 19t 2 s 4 + 48t 3 s 3 + 4t 4 s 2 — 24t 5 s — 81 6 

N c 

2 s 4 t(s + 1) 

4 C F 

(4) 

(s + 2 1) 2 

N c 

s 2 

2 C F 

(5) 

(1 + 2 1) (2s 3 - 3 ts 2 - 2 t 2 s + 2f 3 ) 

N c 

2 sH 

4 C F 

(6) 

(s + 2 1) (s 3 — 7 ts 2 — 8 Ps — 2P) 

N c 

2 s 3 (s + 1) 

4 C F 


Table 2: Expressions for the gg —> gg subprocess corresponding to diagrams (l)-(6) of Fig. [5j 
hence containing only 3-gluon vertices, in gauge (4.39) with non-vanishing 4-gluon vertex con¬ 
tributions. 


In order to add the 4-gluon vertex contribution and obtain a full result for the Di coefficients, 
let us consider a general 4-gluon amplitude, shown on the left hand side of Fig. [7j A 3-gluon 
vertex brings a single SU(N) structure constant factor. Each amplitude in Fig. [5] consists of 
two 3-gluon vertices and that results in three possible color factor products 


c — jaica4yca2d3 ^ _ jd\a.2C c — jaxa^c jca^a2 


(4.40) 


for the amplitudes with a gluon exchange in the t-. s- and u-channels, respectively. Each of the 
above amplitudes can now be written as 




(4.41) 


where i is either t. s or u, Ci is a color factor from Eq. (4.40), and is a corresponding 
kinematic expression. The 3 g superscript means that only 3-gluon vertices are involved in the 
given amplitude. Similarly, for the conjugate amplitudes, following the notation of Fig. [7J we 
have 

- _ fbicbn yeb2&3 ^ — ybi&2Cyc&3&4 - _ jbib^c jcb4,b2 ^ 


That allows us to identify the color coefficients of the 3-gluon diagrams of Fig. [5] and write them 
in a compact form 


(1) <H- c t c t , (2) c u c u , (3) c t c u , 

(4) c s c s , (5) <->■ c s c t , (6) <->■ c s c u . 


(4.43) 


The contributions from diagrams with 4-gluon vertex are depicted in Fig. [6| where the 

first row shows the 4-gluon vertex amplitude squared, and the second row gives the interference 
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A 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 


N c (2s 4 + 2 s 3 t + 3 s 2 P + 8si 3 + 6f 4 ) 
C F s 2 f 2 

1V C (s 4 + 4s 3 i + 15 s 2 t 2 + 16st 3 + 6f 4 ) 
C F s 2 (s + t) 2 

N c (s 4 + s 3 t + 7s 2 t 2 + 12 si 3 + 6f 4 ) 
2 C F s 2 i(s + f) 

A c (s + 2t) 2 

Cps 2 

N c (s — 2t)(l + i)(3 + 2t) 

2 Cps 2 t 

N c t(s + 2t)(3s + 2t) 

2 Cps 2 (s + t) 


Table 3: Full expressions for the diagrams including three-gluon and four-gluon vertex contri¬ 


butions in the gauge (4.39). 


terms with the three types of M 39 amplitudes from Eq. (4.41). A 4-gluon vertex amplitude 


contains all three color factor products of Eq. (4.40) at once 


M 49 = c t Af 9 + c s A 49 + c u Ai 9 . 


(4.44) 


Therefore, all the contributions from Fig. [6] can be represented in the basis of the color factors 
defined in Eq. (4.43). This allows us to distribute all the pieces of diagrams from Fig. [6] over the 
six A expressions, needed to calculate the hard factors (4.33)-(4.35), according to their color 
factors. Hence, the full expressions are 


A 

A 

A 

A 

A 

A 


A (hf + 2A 4 t 9 A 3 t 9 + A?A?) , 

C 2 (h^ + 2A 4 U 9 A 3 U 9 + A 4 U 9 A 4 U 9 ) , 

a ( 4 3) + 4 9 -K 9 +A 9 A 3 u 9 +-K 9 A% 9 ) , 

C 4 (4 3) + 2 A 4 s 9 A 39 + A'J'A'J 1 ) , 

C 5 (4 3) + At 9 A 4 / + At 9 A 3 s 9 + A 4 s 9 A 3g ) , 
C 6 (4 3) + A 49 At 9 + At 9 A 39 + Aj 9 A 3 u 9 ) . 


(4.45) 

(4.46) 

(4.47) 

(4.48) 

(4.49) 

(4.50) 


The results for As in the gauge (4.39) are summarized in in Table[3j Plugging those expressions 
into the hard factor definitions (4.33)-(4.35) leads to the results identical to Eqs. (4.36)-(4.38). 


We have already seen that not all o f the six hard factors that arise in the gg —>• gg subpro¬ 
cess are independent. As shown in Eq. (|4.35 ), the expressions for H gg -+g g , H g 4 ^ gg , an d 


H gg -^gg differ only by numerical factors. On top of that, when examining further Eqs. (4.33), 
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4.34) and (4.35), we see that the hard factors Hgg ^ gg , Hg 2 g^ gg and H g e g L. Sff are linearly depen¬ 


dent, that is 


h ( 6) = H (i) , h (2) 

99^99 99^99 99^99 ' 


(4.51) 


Hence, the cross section for two-gluon production from Eq. (4.29) can be written in a much 


simpler, factorized form, with only two hard factors and two gluon distributions 

„,2 


da pA ~^" x ^ 

d 2 P t d 2 k t dyidy 2 (aq x 2 s) 2 Xl ^9/p( Xl ^ ) 


3>R) + #( 2 ) /y( 2 ) 

99~*99 99~>99 99^99 99^99 


(4.52) 


In this channel, the new gluon TMDs, <]?gg ^ gg , are defined as the following linear combinations 

r x*(l) x*( 2 ) 77 ( 6 ). 

ui j~gg , Jgg , • • • , *rgg • 


<S(!) _ 1 Tr(l)_ Tt(3) 1 _tt(4) _i_tt(5) 1 tt(6) 

99^99 2 \ 99 N 2 99 N 2 99 N 2 99 99 I ’ 


<p( 2 ) 

99^99 


= 7t(2) _ ^ tt(3) 1_L_ tt(4) 1_ - p { 5) 1 tt(6) 

99 jy 2 99 ~ jy2 J 99 ~ jy2 J 99 T J 99 ’ 


and the new hard factors are: 


fti 1 ) = 2 and K^ = H^ 

99->99 99-^99 ’ u 99^99 99^99 ' 


(4.53) 

(4.54) 

(4.55) 


The explicit expressions are given in Table [T] We note, that the above simplification occurs 
naturally when utilizing gauge invariance from the start, as we will show in section [6j 

Finally, we point out that, in the large-A c limit, all the distributions that were introduced 
in this section, F^ F gg \ F g l\ F g ^\ and F gg \ can be written in terms of xG ^ and xG^ 2 \ and 


equivalence of formulas (4.13), (4.25) and (4.52) with CGC results is obtained 113 


Let use conclude that this part of our work brings two improvements to the current state of 
the art for the TMD factorization in forward dijet production. First of all, we have obtained 
finite-A 7 C corrections to the hard factors of Ref. 13 . More importantly, however, we have 


eliminated the redundancy in the number of gluon distributions needed to write a factorization 
formula for this process, which now takes the compact form 


a Q 


dg-pA— >-dijets+X 


d 2 P t d 2 k t dyidy 2 (aq x 2 s) 2 


fa/ P ( x U» 2 )Yl ^ag^fcdP 


1 


a,c,d 


i= 1 


s (0 _ 

ag^cd 1 + 5cd 


(4.56) 


with only two gluon distributions and two hard factors required in each channel. Note that, as 


we shall discuss now, the incoming, small-X gluon is kept on-shell. Eqs. (4.56) will be further 
generalized to the case of the off-shell gluon in Section [5] 


4.4 The \k t \ Q s limit 

Finally, let us consider the limit \kf \ Q s . This is the dilute limit considered in Section [3j with 
the extra requirement that \kt\ -C |T)|, needed for the validity of those formula. In that limit, 
the transverse separation between the field operators in the definition of the gluon distribution 
is restricted to values much smaller than the distance over which the Fourier integrand varies, 
and the £ dependence of the gauge links can be neglected. As a result, they simplify, and all 
the Fa g distributions coincide, except F^ g which vanishes. In terms of the ^g^lcd functions, 
all six distributions also reduce to that one gluon distribution, which can therefore be identified 
with F g/A /7 T. 
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Then, for all channels, one can easily sum the surviving hard factors. In terms of diagrams, 
we always obtain D\ + D 2 + 2 D 3 + D\ + 2 D 5 + 2Dq, meaning that we recover the collinear matrix 
elements. Indeed we have (noting that Hggl^gg + Hgg\g g + Hg^\ gg = 0): 


#(!) _l #(2 ) _ 1) I ^'(2) 

99^99 ~ 99^99 99^99 ~ 99^99 


_ + - + K ^ - 

99^99 ^ 99^99 99^99 ~ 99^99 


+ H (6) = K^> + A' (2) 

99^99 ^ 99^99 99^99 ^ 99^99 


S 2 + U 2 


Cf S 2 + U 2 1 


N r 


t 2 

1 t 2 +u 2 
2N C tu 
2N C (s 2 — tu) c 
C F s 2 i 2 U 2 


= -zl>f 


sn g 
1 P + u 2 
2C f s 2 


99 ^ 99 1 i 


(4.57) 


1 


= 7zl-M, 


99^99 


— ^IM-gg-tqql ,(4.58) 
(4.59) 


sr 

12 


Therefore, we recover the HEF formula (2.9), except that, due to the \kt\ -C \Pt\ limit, the 
matrix elements are on-shell: the transverse momentum of the incoming gluon, kt , survives only 
in J~ g /A■ I n other words, we recover the standard high-|T)| limit: 


^Q-pA-^ dijets+X 

dy\ dy 2 dP?dk 2 


£r 

a.c.d 


+ dcd 


x lfa/p{ x i i k 


da. 


ag^cd 

dt 


Fg/A{x 2 , k t ) , 


(4.60) 


with daag^cjjdt = \M. a g^cd\ /[^’ k (xiX 2 s) 2 ] 1 and where P g /^(x 2 , h) can be identified with 
d/dk 2 X 2 f g /A{ x 2 ,k 2 ), the derivative of the integrated gluon distribution. 

In the following section, we shall restore the kt dependence of the hard factors. This will 
extend our formulas such that they recover the full HEF formula when the dilute limit is con¬ 
sidered. As a result, we will obtain a unified description, valid for generic forward dijet system 
with \pit\, \p 2 t | Qs > without any additional requirement on the magnitude of the transverse 
momentum imbalance kt■ 


Unified description of forward dijets in p+A collisions: 
TMD factorization with off-shell hard factors 


We shall now generalize the hard factors that enter the TMD factorization formula ( |2.10[ ) to 
the case with one of the incoming gluons being off the mass shell, as illustrated in Fig. [ 8 ] As 
it has been already stated, the motivation to include the offshellness is to be able to allow for 
configurations where the dijets are produced at any azimuthal angle (of course before application 
of a jet algorithm that will suppress very small angles and hence render the results finite). 

As can be seen in Fig. [9] (as an example we chose only purely gluonic matrix element but the 
same structure occurs for the other channels), the on-shell matrix element misses substantial 
contributions when the jets are produced at small angles near A<p = 0 and at small rapidity 
differences AY = \y\ — 2/2 1 — 0. In such configurations, the matrix element develops a structure 
that is divergent and it is suppressed only by a jet algorithm, which has to be applied in order 
to ensure two-jet configurations 16 . The matrix elements squared we are after, i.e. gg* -A gg, 
qq and qg* —> qg, can be extracted from the high energy limit (or eikonal limit) of 
qgg and qg —>• qqq and qq 1 -A qq' g 35 . In this approach the quark q is an auxiliary 


99 
99 

line to which the initial state off-shell gluon g* couples eikonally. 

The high energy factorization is a direct procedure where one uses the standard Feynman 
rules for all vertices and color factors, and fixes the light-cone gauge for the on-shell gluons, 
using a gauge vector given by the longitudinal component of the off-shell, initial-state gluon’s 
momentum. In particular, if we apply the high energy factorization to the process we are after, 
we set the gauge vector to n = pa, where pa is the target four-momentum, as defined in Fig. |T] 
and Eq. (2.2). Furthermore, the prescription is to associate with the off-shell gluon a longitudinal 
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Figure 8: Four-parton amplitude with the incoming, small-a;, off-shell gluon, 
polarization vector, called nonsense polarization |l], of the form 2 


6 ° = 
H 1 


iV2 


X2 


I h 


~PA[i ■ 


(5.1) 


As elaborated in Ref. [10 [, longitudinally polarized gluons provide the dominant contribution to 
the cross section in the high energy limit. In the square amplitude, this leads to the polarization 
tensor of the form 110 


e°e°* 


-2xl 


Pa /.i Pa v 


(5.2) 


In the above, X 2 = k^ 1 / pauP u , which follows directly from the definition in Eq. (2.6). The sum 
over polarizations of the on-shell gluons takes the standard form, with the gauge vector given 
by PA 

_ PA^iQu + %PAv 

2.V, -ft- q , np > 


(5.3) 


where, depending on the channel, q = p, Pi or p 2 , c.f. Eq. (4.39). 


Let us note that the procedure outlined above defines the hard process in a gauge invariant 
manner only when a special choice for polarization vectors of the on-shell gluons is taken. In an 
arbitrary gauge, for internal and external gluon lines, more sophisticated methods have to be 



longitudinal component of the off-shell gluon 

S = ( X2PA+P ) 2 = 

t = ( X 2 PA ~ Pi) 2 = 
u = ( X2PA~P2) 2 = 


I Pt \ 2 


z(l-z) 
—zs , 

-(1 - z)s, 


+ \k t \ 2 = Xi X 2 s, 


which are related via the equation 


s T t T u = 0 . 


(5.4a) 

(5.4b) 

(5.4c) 

(5.5) 


4 The \/2 factor in Eq. (5.1) follows from a convention. It allows for use of the on-shell-like factor ^ in averaging 
over polarization, while calculating matrix elements squared, even in the case of the off-shell gluon, where the 
actual number of polarizations in the high energy limit is 1. 
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Figure 9: Matrix elements squared for gg —> gg scattering with p t i = pt 2 = 4 GeV and a s = 0.2. 
Left: the on-shell case. Right: the off-shell case. AY and Acj) are, respectively, the differences 
in rapidity and azimuthal angle of the two outgoing gluons. 


In the on-shell limit, —> 0, the variables defined above recover the standard Mandelstam 

variables from Eq. flTTl ) 


lim (s — s) = 0 . 
|fe|->o 


lim (t — t) = 0 , 
\kt\-tf) 


lim (u — u) = 0 . 
ILKo 


(5.6) 


As a consistency check, we have verified that, for all three subprocesses, the off-shell am¬ 
plitudes that shall be used to build the hard factors in the remaining part of this section are 
identical to those first calculated in Ref. [ll . 

From this point onwards, we shall discuss our results only in terms of the new K® hard 
factors and the new factorization formulas from Eqs. ( 4. 13| ) , ( |4.25[ ) and ( |4.52 ). The results for 
the old hard factors, H^\ in the off-shell case are given in Appendix [a| for completeness. 


5.1 The qg* —> qg channel 


The off-shell hard factors for this channel are obtained using definitions given in Eq. ( |4.10[ ) and 
then Eqs. (4.6) and (4.7). The corresponding Dj expressions are collected in Appendix [A] in 
Table [H The two hard factors read 


K. 


(i) 


s 2 


u 2 


Q9 ~*Q9 


2ttsu 


UU 


SS — tt 

N? 


K, 


( 2 ) 


Cf s (s 2 + n 2 ) 


99 —*99 


N r 


ttu 


(5.7) 

(5.8) 


In the limit \kt\ —> 0, simplification given by Eq. (5.6) occurs and the above formulas mani¬ 
festly recover the on-shell results from Table [l] 


5.2 The gg* —> qq channel 


The off-shell hard factors are obtained using definitions given in Eq. ( |4.28| ) and then Eqs. (4.18), 
(4.19) and (4.20). The corresponding Dj expressions are collected in Appendix |A] in Table [9j 
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The two hard factors take the following compact form 


K 


K. 


(1) 

gg*^qq 

( 2 ) 

99 * ^qq 


1 t 2 + u 2 
2 N c sstu 

1 


t 2 + u 2 


4N%Cf sstu 


[uu + tf] , 

[uu + tt— ss] 


(5.9) 

(5.10) 


Again, following Eq. (5.6), it is manifest that the above hard factors reduce to those given in 
Table [lj in the limit \kt\ —> 0. 

5.3 The gg* —> gg channel 

In the gauge chosen for our calculation, all the squared diagrams and interference terms that 
involve a 4-gluon vertex are identically zero. The corresponding DjS are given in Table 10 
of Appendix |A} Using the combinations from Eqs. (4.33)-(4.35) and then the definition from 
Eq. (4.55) leads to the following set of the off-shell hard factors 


K. 


(i) 


2N C (s 2 — tu) 2 


K. 


99 —*99 
( 2 ) 


99 ~*99 


Cf ttuuss 
N c (s 2 — tu) 2 
Cf ttuuss 


[uu + tt\ , 

[uu + it — ss] . 


(5.11) 

(5.12) 


The on-shell limit is again manifest, with the above equations reducing to those from Table [T] 
as \kt\ —> 0. 

6 Helicity method for TMD amplitudes 

(z) 

In the preceding sections, the hard factors accompanying the gluon densities Tag were calcu¬ 
lated from the squared diagrams presented in Figs. [3](6j This procedure has certain drawbacks, 
especially when one would like to consider more complicated processes. For multiparticle pro¬ 


cesses, 


the color decompositions and helicity method 26, 49 are now considered as the most 


effective ways to deal with them. Moreover, it is not obvious how the gauge invariance comes 
into play for the separate diagrams from Figs. [3]{6] contributing to the hard factors. In the color 
decomposition method, the so-called color ordered amplitudes are gauge invariant from the start 
and one can use them directly to construct hard factors. 

In view of the above, and to cross-check the results from Section[5j we will give an alternative 
procedure to obtain the factorization formulas with off-shell gluon. To this end, we shall need 
TMD gluon densities corresponding to color decomposition of amplitudes and the color-ordered 
amplitudes themselves. 


6.1 Color decompositions 


Let us recall some basic facts about the color decompositions. We refer to 26,49 for more 
details. 


We first consider a gluon amplitude A4 ai- " ajv ,..., ^, where ai,..., ajy are the 


ex¬ 


ternal, adjoint color quantum numbers, the is a polarization vector for a gluon i having 
momentum ki and helicity Aj = ±. The fundamental color decomposition reads 


j^a i...a/v / gAl 


1 ,...,4")= £ Tr(t a >f 

CtESjv — 1 


°2 


. t a ”N) M (l Al , a. 


•^cr2 —.^crN 

2 * * * 5 a N 


( 6 . 1 ) 
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where the sum is over a set SV-i of all non-cyclic permutations of {1,, N}. The coefficients 
of the expansion define color ordered - or dual - amplitudes. They possess several useful 
properties. First of all, they are gauge invariant. Second, there are certain relations between dual 
amplitudes. Indeed, the following adjoint color decomposition involves only (IV — 2)! different 


amplitudes 50 


M 




> °N 


= £ {F a ^...F a °^) aiaN M(l x \o } 2 


(T£Sjj —2 


,CT 


MJV - 
N—l 


,N* 


( 6 . 2 ) 

where ( F a ) bc = f abc . _ 

Consider now an amplitude involving a quark anti-quark pair J\y[ D i a 2 -~ a N-iD N w i iere 
Dj are the color and the anti-color of the quark and the anti-quark, respectively. The color 
decomposition reads 


M D ia2 ...a N _ lDN (x l ,4\...,£^ r S 1 1 ,X N ^ = 

E • ■ ■ t a " N ~ l ) N M (l A| , ..., , AT a ") ■ (6.3) 

Now Ai and Atv are helicities of the quark and the anti-quark. For amplitudes involving more 
quark anti-quark pairs the decomposition is more complicated and we refer to [261 for details. 

It is important to note that the above color decompositions work also for the case when one 
of the gluons is off-shell. 


6.2 Gluon TMDs for color ordered amplitudes 

Let us now fold the gluon TMDs corresponding to the color ordered amplitudes squared, as de¬ 
fined in the previous subsection. We constraint ourselves to the 2 —> 2 processes case considered 
in this paper. 

Let us first consider the g (£q) g* (k\) —> g (k^) g (k 2) process. For the purpose of this and 
next subsections we have assigned a new set of momenta to the partons. This assignment differs 
from the one used before but it is more convenient when dealing with color ordered amplitudes. 
The correspondence is achieved by the following relations: k± «-)• k, &2 pi, k^, <->• p- 2 , &4 -H- p. 
Moreover, for the off-shell momentum we adopt a notation 

ki = n\ + kx ■ (6-4) 

The color decomposition of the four gluon amplitude reads 


\ja1a2a3a4 f \ 2 f \3 p A 4 \ _ f f kA 

J l gg*^gg ) fc 3 > fc 4 J ~ Ja 1 a2Cjca3a4 


1 * 9 A 2 

99*^99 \ 1 J z 


,3 As ,4 A4 ) 


+ f a\ci3cfca 2 cL4 ■M-gg*^gg I 4 


3 Aa , 2 A 2,4 A 4^ , (6.5) 


where ni is placed for the off-shell gluon instead of a polarization vector (in fact it plays a 
similar role). As far as dual amplitudes are concerned, we indicate the off-shell gluon by a star. 
In Table [4j we calculate the gluon TMDs that correspond to the color structures exposed in 
(6.5) (after squaring). They agree with the gluon TMDs calculated in [l2| and listed in rows 1 
and 3 of Table 8 of 1121. That table defines one more gluon TMD (the row 2) which however 
is redundant. Clearly, the color decomposition (6.5) gives all the necessary color structures and 


already incorporates the gauge invariance. In summary, the two gluon TMD listed in Table [4] 
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color-ordered amplitude squared 

gluon TMD 


Mgg^gg (l* , 2 A2 ,3 A3 ,4 A4 ) | 2 

Mgg^gg (l* , 3 A3 ,2 A2 ,4 A4 ) f ' 

*$-+99= 2Nf(^S ) -2^ ) 

+4t ) +^9+N^S ) ) 

M gg ^ gg (l*,2 A2 ,3 A3 ,4 A4 ) M* gg ^ gg (1*,3 A3 ,2 A2 ,4 A4 ) 
M* gg ^ gg (1*,2 A2 ,3 A3 ,4 ^)M gg *^ gg (D,3 A3 ,2 A2 ,4 A4 ) 

^ 99 = W c ( N c4f- 2-^i? 

+jjg > +4g > +N^jjf) 


Table 4: Gluon TMDs accompanying the color-ordered amplitudes for gg* —> q q pro cess. 


H) ■ ■ 

has been assumed that TMDs are real. The J~ gg distributions are defined in Eqs. (4.15), (4.16), 
(4171) and in Eqs. Ool), (|43lj), Oil). 


color-ordered amplitude squared 

gluon TMD 

\M gg ^ qg (2 A2 ,D,4 A4 ,3 A3 )| 2 
\M gg ^ qg (2 A2 ,4 A4 ,1*,3 A3 )| 2 

aM) —if aj2 jr(l) -zr( 3 )A 

^gg^qq ~~ n£-i y^c^gg Jgg J 

Mgg^gg 2^ , 1* , 4^ , 3^ ) M* gg ^ ^ , 4^ , 1* , 3^) 
M* g ^ q - q (2 A2 , 1* , 4 A4 ,3 A3 ) Mgg^gg (2^ , 4^ , 1* , 3^) 

(p( 2 ) _ _7V 2 jr( 2 ) i tt( 3 ) 

^99^99 ~ JV c- r 99 + '99 


Table 5: Gluon TMDs accompanying the color-ordered amplitudes for gg —> qq process. 

(i) r 

been assumed that correlators are real. The T gg distributions are defined in Eqs. (4.15), 
and (|4.ITh . 


t has 

416l 


are the only relevant TMDs and correspond to the two independent gauge invariant amplitudes 
squared and their interference. 

Now, let us turn to the g (Aq) g* ( k\ ) —»• q (fc 3 ) q (A' 2 ) process. The color decomposition reads 

+ 2 i *,4 A4 ,l*,3 A ») . (6.6) 

The gluon TMDs corresponding to the color structures appearing after squaring this equation 
are gathered in Table [5j They correspond to rows 1 and 5 of Table 7 in |l2j. Again, we have 
only two independent TMDs that are needed. 

For the process q (& 4 ) g* (k\) —> q (/c 3 ) g (fo), the color decomposition reads 

“f 4 (A 3 ,n 1 ,4 2 ,A 4 ) = (t a 't a2 )o 3 D A M qg ^ qg (s A3 ,1*, 2 A2 ,4 A4 ) + 

("% 3 d 4 ^V-^ (3 As , 2 Az , D,4 A4 ) . (6.7) 

For anti-quarks we need to exchange the indices 3 4. The TMDs corresponding to those 

processes are given in Table [6j In general, the TMDs for a sub-process with anti-quarks are 
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color-ordered amplitude squared 

gluon TMD 

Mgg^gg (3 ^ , 1 * , 2 ^ , 4^ ) M* qg ^ qg (3^,2^, 1*,4^) 
M* qg ^ qg (3 A 3 ,1*,2 A 2 ,A 4 ) Mgg^gg (3^,2^, 1*,4^) 
|AV^ 9 (3 A 3 ,2 a M*,4 A4 )| 2 

X 

II 

cn 

Cr 

K s *^ 9 (3 A 3 ,1*,2 A 2 ,4 A4 )| 2 

a>( 2 ) _ 1 ( ttO-) I t \ t 2 t ( 2 )\ 

~ n'*—i \ ^19 ' %r qg J 

Table 6 : Gluon TMDs accompanying the color-ordered amplitudes for qq* — > qq process. I 


(i) . . . 

has been assumed that correlators are real. The J- gg distributions are defined in Eqs. (4.4) and 


Eqs. (4.5). 


different than for quarks, but they turn out to be the same assuming that the correlators are 
real. Again, we end up with only two independent TMDs. 


6.3 Off-shell color-ordered helicity amplitudes 


In Section [5j we have calculated the off-shell hard factors in a specific axial gauge, with pa 
chosen as the gauge vector, and using the high energy projector (5.1). As shown in Ref. 1101, 


such a procedure yields results which are gauge invariant within a subclass of axial gauges with 
the gauge vector n ^ = ap g + where a and b are arbitrary complex numbers. There are also 
methods to calculate gauge invariant off-shell amplitudes in any gauge and choice of polarization 
vectors 134,35,47,48]. In what follows, we shall use those methods and specifically the results 
of 35,481. 


Consider first the gluon amplitudes. For the purpose of this section only we assume all 
momenta to be outgoing. For the non-vanishing helicity configurations, in the helicity basis, we 
have 


Mr™(i*>2-,3 + ,4 + ) 
■Mg*g^gg (l ,2 + ,3 ,4 + ) 
■Mg*g^gg (l ,2 + ,3 + ,4 ) 


2g Pl (1*2)(23)(34)(41*) ’ 

2 2 (1*3 ) 4 

9 Pl (1*2)(23)(34)(41*) ’ 

O 2 (1*4 ) 4 

9 Pl (1*2)(23)(34)(41*) ’ 


( 6 . 8 ) 

(6.9) 

( 6 . 10 ) 


where we adopted a shorthand notation for the spinor products (ij) = (ki — \kj+) with |fc«±) = 
2 (1 ± 75 ) u (ki), and where pi is a, for our purposes irrelevant, phase factor (see details e.g. 
in 14 8]). We also defined (l*i) = ( n\i ) with n\ being the longitudinal component of Aq, c.f. 
Eq. (6.4). The other remaining helicity configurations can be obtained from Eqs. ( 6 . 8 )-(6.10) 
using CP invariance 


Mgg^gg 3~, 4~) =M* gg .^ gg {l*, 2~, 3 + , 4 + ), ( 6 . 11 ) 

and so on. For the other color ordered amplitude, M. gg *^ gg (1*, 3, 2,4), we need to exchange 
2 -H- 3 in the denominators. 

The above helicity amplitudes can be efficiently evaluated and squared numerically, however 
for the purpose of this paper we shall need analytic expressions. To this end let us introduce 
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[ij] = (ki + \kj —), which, up to an unimportant phase, is a complex conjugate of (ij). Moreover, 
we have the following relation 

(ij)\ji] = {h + kj ) 2 = Sij. ( 6 . 12 ) 

For the products involving m we use the notation 

<!**)[*!*] = (m + ki) 2 = Si*i- (6.13) 


With this, we get for the required amplitudes squared summed and averaged over helicities 

~4 . ~4 , ~4 

4 S lt2 + 3 + § 4*4 


lM gg *^ gg (l*,2,3,4)l =8g 
jM gg ^ gg (l*,3,2,4)j 2 = 8g 


Sl*2$23S34S41* 

~4 i ~4 i ~4 

4 -S pa + gj*3 + s l*4 

Sl*3S32S24S41* 

~4 , ~4 i ~4 

4 S lt2 + 54*3 + 54,4 


M gg *^ ga (1*, 2,3,4) M gg *^ gg (1*, 3,2,4) - 85 (p^g^pg^g^ 

7W* _ (1*, 2, 3,4) (1*, 3, 2,4) = -8 g A ^* 2 + + S '* 4 


[l*2][34](l*3)(24)s 23 S4i* 


(6.14) 

(6.15) 

(6.16) 
(6.17) 


where we have used overlines to indicate helicity summations. The last two interference terms 
enter the cross section as a sum. Therefore, we may simplify it as 


M gg *^ gg (l*,2,3,4)M 


9g ^ gg (1*, 3,2,4) + M gg ^ gg ( 1 *, 2 ,3,4) M gg ^ gg ( 1 *, 3, 2 ,4) 

4 (®1*2 + sf*3 + ^ 4 * 4 )(§24^1*3 — S 2 3Si*4 + S34Si* 2 ) 


= — 8 y 


Si* 2 S34Si* 3 S24-S 2 3'S41* 


(6.18) 


where we have used 


[1*2][34](1*3)(24) + (1*2)(34)[1*3][24] = (m- |^ 4 # 2 |ni-) + (m- , (6.19) 

and applied = s tJ — jf-j'jji, a few times. The amplitudes for the on-shell limit are simply 
obtained by dropping the star in 1 * so that the spinor and the scalar products will be with k\ 
instead of n±. 

Now let us turn to processes with quarks. We will give only amplitudes for g (£q) g* (&q) —> 
q (A 3 ) q (k 2 ) process, as all the other can be obtained by the crossing symmetry (taking care of 
the proper color flow when crossing). We have 

A4,^(3-,l*,4+,2+)=2(6.20) 

< 3+ . i*,4 + ,2~) = 2g 2 * <6 ’ 21> 

We note that the above formulas have never been published in the literature and are given here 
for the first time. 

Similar as before, the two remaining helicity configurations can be obtained thanks to CP 
symmetry. For the color ordered amplitudes with 1 and 4 interchanged, we need to make a 
replacement 1 -H- 4 in the denominators. The amplitudes squared and summed over helicities 
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read (the helicity averaging factor is included) 


•^99*->99 
•^99*->99 


|-^ 99*->99 (3) 1 *) 4, 2 ) | 2 = 
|M gM (3,4,l*,2 )| 2 = 
(3, 1 *, 4, 2 ) M* gg *->qq (3, 4, 1 *, 2 ) 
(3,1*, 4, 2) Mgg*^. q q (3,4,1*,2) 


V £i^k+iM_ 

Si*4S34S23 

V^i«, 

Si*4S24S23 

4 Sl* 2 Sl*3 (si*2 + ®1* 3 ) 

" ~~ 9 (21*)(43)[31*][42]S 23 S 41 * ’ 

_ _2 4 Sl*2Sl*3 (^4*2 + ^1*3) 

9 [21*][43](31*)(42)s 23 S4i*' 


The sum of the last two interference terms simplifies to 


( 6 . 22 ) 

(6.23) 

(6.24) 

(6.25) 


-^ 99*->99 (3, 1 , 4, 2) ■M-gg*^qq (3, 4, 1 , 2) + A4 gg *^ qq (3, 1 , 4, 2) M gg *^ qq (3, 4, 1 ,2) 

4 Sl*2Sl*3(sf» 2 + ^1*3)(®24'®1*3 — S23«l*4 + 53451*2) 


= -2<f 


Sl*2'S34-Sl*3'S24S23-S41* 


(6.26) 


In order to obtain amplitudes for q (£ 4 ) 5 * (&i) —)> q (£ 3 ) g (^ 2 ) we can use the crossing symme¬ 
try. Specifically, we can obtain \Mqg*-+qg (3,1*,2,4)| , \M. q g*^ fqg (3,2, l*,4)j and interference 


terms by making replacement 2f>4in Eqs. ( |6.23[ ), ( | 6 .22 [ ) , ( |6.26[ ) respectively. 


6.4 Hard factors from color-ordered amplitudes 

Having computed the color ordered amplitudes it is now straightforward to calculate the hard 
factors J\(d. Let us note, that it is the K W hard factors that appear naturally within the 
color-ordered formalism, not the factors. It also comes naturally that there are two hard 
factors and two TMDs per each channel, so the the factorization formulas can be written in a 
unified form: 


dijets+X 

d 2 P t d 2 ktdyidy 2 


a c 


(xiX2S) 2 ^ 
a,c,d 


•^1 fa/pi^l » M ) £ ^ag*^cdP 1 


i=l 


5 (0 1 

ag^cd 1+ ^ cd 


(6.27) 


where a,c,d are the contributing partons. The explicit expressions for the generalized gluon 


TMDs 


(0 

ag^cd 


are listed in Tables 


6 The hard factors K 1 were already given in Section 5 (we 
collect them in Table [7] for convenience). In the context of this section, they are obtained by 
multiplying the left column of Tables [4j{6] by the corresponding color factors and combining the 
cells that belong to the same generalized TMD. More precisely, we have 
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K {1 1 
99 ->99 


9 


4 


K^l 

99 *->99 


1 N*C F 

(2N c C f ) 2 2 
1 N^Cp 

(2 N c C f ) 2 4 


(| -^ 99*—>99 (1 >2,3, 4) | + \j\Agg*^gg (1 , 3, 2, 4) | ^ , 
(Mgg^gg (1*, 2, 3, 4) M* gg .^ gg (1* , 3, 2, 4) + C.C.) , 


(6.28) 

(6.29) 


for pure gluon channel, and 


9 
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K (1 l - = 

99 ->99 
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4 


k (2 1 - = 

99 “>99 


1 

(2 N c C f ) 2 
1 

(2 N c C f ) 2 


N c Cp 

-C F 

2 


(|AV ^99 (3,1*,4, 2)| 2 + | M gg ^ qq (3,4, 1 *, 2 )| 2 ) 
(Mgg^qqi 3,1*, 4,2) A4 gg *^. qq (3,4,1*,2)+ c.c.) , 


(6.30) 

(6.31) 
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Table 7: The hard factors accompanying the gluon TMDs <& 


(0 

ag^cd’ 


for gg* —> qq channel. For the qg* —> qg sub-process we need to use the crossing symmetry as 
described in the preceding section. We have 


o 4 

y ^gg* 


qg*^qg 


( I -Mgg^qq (3, 1 ,4, 2)| 


2 CpN* 
Ci 
2 


2^4 


r> f ■Mgg*^,qq (3, 1 , 4, 2) Aigg*^ q q (3, 4, 1 ,2) C.C.^^ | , 


= 5^W^(-|M 99 .^ 5 (3,4.r,2)| 2 ) 


2f>4 


(6.32) 

(6.33) 


In all the formulas above, the first color factor comes from color averaging. The minus signs in 
front of the amplitudes in (6.32), (6.33) come from the crossing of a fermion line. Table [7] is 
easily recovered using the following relations of §ij to the kinematic variables from Section 


S23 = Sl4 = S, S34 = S12 = t, S24 = S13 = U , 
Sl*4 = s, Sl*2 = t, Sl*3 = U . 


(6.34) 

(6.35) 


7 Conclusions and outlook 


Dijet production is one of the key processes studied at the LHC. Requiring the two jets to be 
produced in the forward direction creates an asymmetric situation, in which one of the incoming 
hadrons is probed at large x, while the other is probed at a very small momentum fraction. This 
kinematic regime poses various challenges, one of the biggest questions being the existence of 
a theoretically-consistent and, at the same time, practically-manageable factorization formula. 
The standard collinear factorization is not applicable in this case as the dependence on the 
transverse momentum of the low-x gluon in the target, kt 1 cannot be neglected. 

In the limit where the jets’ transverse momenta \pit\, \p 2 t\ \h\ ~ Qs, with the latter being 

the saturation scale of the target, an effective transverse-momentum-dependent factorization 
formula for forward dijet production has been derived in Refs. |13,l4j and it has been shown 
to be consistent with the CGC framework. On the other side, the high energy factorization 


approach 10 has been also successfully applied for studying forward dijet production at the 
LHC. In this paper, we have examined the theoretical status of the HEF approach in the context 
of forward dijet production at hadron colliders and reconciled it with the TMD factorization by 
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creating a unified framework valid in the limit \pu\-, \P 2 t\ ^ Qs with an arbitrary value of \kt\, as 
long as it is allowed by phase space constraints. In particular, we have shown in Section [3] that 
the HEF formula is indeed justified in the kinematic window of \pu\, \p 2 t\ ~ \kt\ 3> Q s , where 
it was explicitly derived from CGC for all 2 —>■ 2 channels. This limit corresponds to the dilute 
target approximation hence no non-linear effects are expected. 

The second major result of our work is an improvement of the effective TMD factorization 
for forward dijet production, first derived in Ref. jl3], by taking into account in Section [4] 
all finite-A^ c corrections, as well as generalizing the factorization formula to the case with an 
off-shell incoming gluon in Sections [5] and [6j In addition, we were able to simplify the TMD 
factorization formula by reducing the number of gluon distributions to two independent TMDs 
for each channel. The main results of this part of our study are summarized in Eq. ( |6.27 ), 
which gives the new TMD factorization formula, as well as in Table [7j where we collect all the 
off-shell hard factors. The corresponding gluon distributions are given in Tables [4j [5] and [6j The 
above results were obtained with two independent techniques: a traditional Feynman diagram 
approach and helicity methods with color ordered amplitudes. The improved TMD factorization 
formula (6.27) encapsulates both the result of Ref. 1131 and the HEF framework as its limiting 
cases. 

The results obtained in this paper open several avenues for future research that we plan 
to follow. First, a natural next steps will be to use Eq. (6.27) for phenomenological studies. 
That shall require some input for the six gluon TMDs which may be difficult in 

a general case. But in the large- N c limit, they can all be written in terms of just two functions: 
xG^ l \x,kt) and xG^ 2 \x,kt), which in turn can be evaluated within certain models, as in |5|. 

Another line of possible extension of our framework is to supplement it with high-|F)| effects 
such as Sudakov logarithms or coherence in the evolution of the gluon density. Essentially, this 
can be done by adding a p 2 dependence to the unintegrated gluon distributions 29 31,51 - 54 


The equations that combine such effects with the small-® evolution (55,56 show a nontrivial 
interplay between the non-linearities and the p 2 dependence and this may, in particular, weaken 
the saturation effects. At the linear level, the so-called single step inclusion of the hard-scale 
effects (as demonstrated in |17|) helps in the description of forward-central dijet data, therefore 
this direction seems to be relevant in order to provide complete predictions. Furthermore, first 
estimates of azimuthal decorrelations of the forward-forward dijets in the HEF framework, with 
inclusion of hard scale effects and non-linearities, show that they are of similar relevance for this 
process 33]. 

Last but not least, it remains to be proved that the large logarithms generated by higher- 
order corrections can indeed be absorbed into evolution equations for the various parton dis¬ 
tributions (and jet fragmentation functions) involved, and potentially for additional soft fac- 
This limitation however is not specific to our work, the same is true at the level of 


57 


tors 

the TMD and HEF regimes independently. In the former case, it is known that TMD factor¬ 
ization generically does not apply for dijet production in hadron-hadron collisions 22,24 . It is 


nevertheless expected that, in dilute-dense collisions, initial state interactions originating from 
a dilute hadron do not interfere with the intrinsic transverse momentum and thus factorization 
may hold, although there is no formal proof of this statement yet. 

In addition, even though it was possible to write formula (4.56) in terms of just two TMDs per 
channel, this simplification may not survive after small-® evolution is included, as, in general, the 
non-linear equations mix the original J-jig functions. For instance, xG^ does not obey a closed 

We 
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equation and, contrary to what happens with xG^ 2 \ the large- N c limit does not help 
note that any equivalent linear combination of the gluon distributions, such as (2.10) and (4.56), 
is equally valid, and it may turn out that some alternative choice allows one to write the 
evolution equations directly in terms of TMDs. By contrast, it is also possible that the inclusion 
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of small-x evolution can only be achieved within the full complexity of the CGC, meaning that 
the Q s ~ \kt.\ -C |-Pf| limit, which allows one to avoid the quadrupole operator in (3.10) and 
express the cross section in terms of gluon distributions, may not help when small-x evolution 
is considered. 

In the HEF regime, the issues are different. The Q s <C \kt\ ~ \Pt\ limit makes things simpler 
from the point of view of small-x evolution, since non-linear effects can be neglected. However, 
the off-shellness of the hard process is not neglected and thus the standard power counting of 
the twist expansion becomes useless. One must then resort to different methods, such as those 
Any progress towards an all-order proof of either HEF or TMD factorization for 


of Ref. 59 


forward dijet production in dilute-dense collisions will naturally carry over to our improved 
TMD factorization formula (6.27) that combines both regimes. In the meantime, our results 
represent a viable alternative to CGC calculations, equivalent to them in the kinematic regime 
appropriate for dijets Q s <C \Pt\ but more practical. 
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A Off-shell expressions 


m 


In this appendix, we gather all expressions corresponding to the D{ diagrams from Fig 
the case where one of the incoming gluons is off-shell. All calculations were preformed in the 
axial gauge discussed at the beginning of Section [5j with the axial vectors for the on-shell gluons 
set according to Eq. (4.39). 

For completeness, we also give here the results for the “old” hard factors defined in Eqs. (4.6), 


( |4.7[ ) ( |4. 18 ), ( |4.19[ ), ( |4.20[ ), ( |4.33[ ), ( |4.34[ ) and ( |4.35[ ), in the case with off-shell incoming gluon. 
Table [8] gives the D{ expressions for the subprocesses qg* —> qg. The two hard factors in this 


channel read 


H, 


(i) 


qg —*qg 


s 2 + u 2 
2 stu 


u — 


tt 

N^u 


(A.l) 


H, 


( 2 ) 

qg*^qg 


S (s 2 + U 2 ) 


2 utt 


(A.2) 


In the limit, \kt\ —>• 0, simplification given by Eq. (5.6) occurs and the above formulas manifestly 


recover the on-shell results from Eqs. (4.8) and (4.9). 


The corresponding Di results for the gg* —> qq subprocess are given in Table [9] The three 
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qg* qg 

A 

(i) 

2 tu + P + 2 u 2 

7 (s + 7 + u) 

(2) 

Cf u (2 tu + P + 2 u 2 ) 

N c ut (s + t + u) 

(3) 

(2 tu + P + 2 u 2 ) (t (s + 7) + u (s + u)) 

4 tut (s + 7 + u) 

(4) 

Cf (t + ^0 {p'tu H - 
N c §t (s + 7 + n) 

(5) 

(27ft + T 2 + 2tZ 2 ) (7 (s — 7) + u (s + «)) 
4s77 (S + 7 + u) 

(6) 

1 (27u -)- P + 2n 2 ) (7 (s + 7) + iZ (s — «)) 

A? 4stt7 (s + 7 + u) 


Table 8: Expressions for the qg* —> qg subprocess with off-shell incoming gluon corresponding 
to diagrams (l)-(6) of Fig. [3] in gauge described in Section [5} 


“old”, off-shell hard factors for this channel take the form 

H m 

H. 


99*^99 

r(2) 


H. 


99 -><?<? 
(3) 


99 —*qq 


1 P + u 2 r A _ 

77 :,- 77 ^- \uu + tt\ , 

4 Cp utss 

1 P + u 2 r_ — A _-, 

—---- \ss — tt — uu\ 

4 Cf utss 

1 P + u 2 

4 N^Cf ut 


(A.3) 

(A.4) 

(A.5) 


Again, following Eq. (5.6), it is manifest that the above hard factors reduce to Eqs. (4.21)-(4.23) 
in the limit \kt\ —> 0. 


Finally, the A expressions for the subprocess gg* —> gg are given in Table 10 and the six 
hard factors read 


^C 1 ) _ 

99*^99 

N c (s 2 — tu) 2 

Cf ttuuss 
N c (s 2 - iu) 2 

rr( 2 ) 

.9.9*793 

Cf ttuuss 
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99*^99 
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uu 


[ft + 1 
[ss — tt — uu ] 


= N 2 h[ 4 } 


= n 2 h[ 5 1 


N c ( s 2 — tu) 2 


(A.6) 

(A.7) 

(A.8) 


The on-shell limit is again manifest, with the above equations reducing to Eqs. (4.36), (4.37) 


and (4.38) as \kt\ -7 0. 
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55* 99 

A: 


(1) 

1 (s + u) (2 su + s 2 + 2 u 2 ) 

N c 2 ts (s + t + u) 


(2) 

1 u (2 su + s 2 + 2 it 2 ) 

N c 2 us (s + t + it) 
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1 (2 su + s 2 + 2n 2 ) (s (s + t) + u (t - 


N%Cf 8tus(s + t + u ) 


(4) 

1 2 su + s 2 + 2 u 2 

Cf 2s (s +1 + uj 


(5) 

1 (2 su + s 2 + 2u 2 ) (s (s — f) — u (t + 

«)) 

Cf 8sis (s + t + u) 


(6) 

1 (2 su + s 2 + 2u 2 ) (s (s + t) +u(t + 

«)) 

Cf 8sus(s + t + u) 



Table 9: Expressions for the gg* —> qq subprocess with off-shell incoming gluon corresponding 
to diagrams (l)-( 6 ) of Fig. [4] in gauge described in Section [5} 
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(1) 

21V C (A + f 2 + u 2 ) 2 

Cf tu(t + u) (s T t T u) 

(2) 

2iV c (A + f 2 + u 2 ) 2 

Cp ut (t + u) (s + t + «) 

(3) 

iV c (tw + f 2 + u 2 ) 2 (f (s + t) + u (s + u)) 

2 Cf tutu (t + u) (s + t + it) 
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2N C ( tu + t 2 + u 2 ) 2 
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N c (tu + f 2 + u 2 ) “ (t (s — t) + ft (s + it)) 
2Cf stfiZ (t + ii) (s + t + u) 
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N c (tu + t 2 + ti 2 )" (t (s + f) + u(s — u)) 

2 Cf suiu (t + u) (s +t + u) 


Table 10: Expressions for the gg* —> gg subprocess with off-shell incoming gluon in gauge 
described in Section [5j The numbering (l)-( 6 ) corresponds to the color structures as defined in 
Eq. (4.43) and each expression contains contributions from diagrams with both 3- and 4-gluon 
vertices. 
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